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Abstract 



> 

^ . 

•/^ , We propose a new form of the frustrated Heisenberg antiferro- 

magnetic Hamiltonian with spatiahy anisotropic exchange parameters 

^^ . J^^ J^^ and Jac extended along the c, a, and a it c lattice directions 

and apply it to describe fascinating physics of copper carbodiimide 
CuNCN in the assumption of resonanting valence bond (RVB) type 
of its phases. These are invoked to explain the intriguing absence of 
magnetic order in CuNCN down to 4 K. We show that the quasiparti- 
cle spectrum of the RVB model of the proposed Hamiltonian has three 
principal regimes: (i) one with two pairs of lines of nodes, (ii) one with 
a pair of lines of nodes (termed as ID- and QID-RVB states), (iii) and 
one with two pseudogaps and four nodal points (2D-RVB). We present 
a complete parameters-temperature phase diagram of the c-a-ca- RVB 
model constructed with use of the high-temperature expansion of the 
free energy. The phase diagram thus obtained contains eight different 
phases whose magnetic behavior includes Curie and Pauli paramag- 
netism (respectively, in disordered and ID- or QID-RVB phases), and 
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gapped (quasi- Arrhenius) paramagnetism (2D-RVB phases). Adding 
magnetostriction and elastic terms to the free energy of the model we 
derive possible structural manifestations of transitions between vari- 
ous RVB phases of the c-a-ca-model of CuNCN. Assuming a sequence 
of transitions between RVB phases to occur in CuNCN while temper- 
ature decreases explains the features observed in the temperature runs 
of the magnetic susceptibility and lattice constants. Confronting these 
with the magnetic susceptibility and strucutre data measured as func- 
tions of temperature in the range between ca. 20 and 200 K we show 
a remarkably good agreement between our theoretical predictions and 
the experiment as reached by ascribing the model parameters values 
which are intuitively acceptable both in terms of their absolute and 
relative magnitudes and of the character of their geometry dependence. 
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1 Introduction 

Recently CuNCN phase had been obtained and a series of measurements had 
been performed of its spatial structure and magnetic susceptibility, electric 
resistivity, heat capacity (all vs. T) [Ij. Although on the basis of analogy 
with other materials of the MNCN series (M = Mn, Fe, Co, Ni) one would 
expect more or less standard antiferromagnetic behavior, it turned out that 
at low temperature this material does not manifest any magnetic neutron 
scattering [2j. The plausible explanation of the latter fact as "absence of local 
momenta" may be, however, misleading since the absence of the magnetic 
scattering means only the absence of the long range magnetic order (LRMO) 
(evanscence of the spin-spin correlation function), not the absence of the 
local momenta themselves. Similar situation can be observed in the crystals 
of Cu-carboxylate dimers where local spins 1/2 do present on each Cu^+ ion, 
but form isolated singlet pairs so that no long range magnetic order does 
exist. On the structural grounds one cannot expect anything like this in 
CuNCN since among the contacts of individual Cu^^ ions one cannot select 
any which would be uniquely strong. This brought usj3[ HI O [6] to the idea 
that the ground state of this material may be related to the RVB state of 
the Cu^"^ local spins 1/2. 

Previously we assumed [31 IH El E] that from the materials structure Fig. [T] 
[IJ one can conjecture an anisotropic triangular antiferromagnetic Heisenberg 
model with the plane corresponding to the ah plane of the structure. This 
was not easy to reconsile with the intuitive picture [7] of the most important 
couplings to be extended in the a and c directions. This brings us to the idea 
to consider more antiferromagnetic couplings. This results in the Heisenberg 
model with the Hamiltonian: 

^5^J.SrSr+. (1) 

r T 

where the translation vector r takes four values Tf, i = 1 -^ 4; Ti = (a, 0);t2 = 
(0,c);r3 = (a, c);r4 = {a,—c) with the interaction of the strength J^ along 
the lattice vector Ti (two neighbors), with a strength J^ along the lattice 
vector T2 (two neighbors as well), and interacion of the strength Jac along 
the lattice vectors t^ and T4, (two neighbors along each). The importance of 
the diagonal Jac couplings in the r^ and T4 directions have been recently re- 
iterated in j8] although it is fairly in a line with the standard considerations 
of [7]. Either interactions along Ti and T2 or those along t^ and T4 taken 





Figure 1: The CuNCN crystal structure and the exchange parameters in- 
cluded in the consideration, (a) shows the overall view on the structure. A 
stronger J^ extends in the c direction; somewhat weaker J^ extends in the 
a-direction. The weakest Jab extends along the 6 ± a directions, (b) shows 
the ac planes of CuNCN. The interactions in all cases are mediated by the 
NCN^~ moieties. Two stronger interactions {Jc and Jac) are mediated by 
the TT-system of NCN^~; somewhat weaker J a is strongly contributed by a 
ferromagnetic cointerpoise terms dependent on the hybridization at the N 
atoms, (c) shows for completeness the ah planes of the CuNCN which are 
not considered in the present paper. 



separately must lead to an antiferromagnetic state. However, when consid- 
ered simulateneously they interfere leading to a frustration not allowing the 
spins to arrange in any LRMO state. For similar systems a variety of RVB 
states have been proposed j9l [10]. Ground state of a similar, but spatially 
isotropic J1J2J3 model have been treated recently by various methods and it 
has been shown that spin-liquid states are very probable [IIj. In the present 
paper we consider in details the RVB states of the above model in the RVB 
mean-field approximation and apply this to analysis of the experimental data 
so far obtained for CuNCN. 

2 RVB mean-field analysis of the model 

2.1 Quasiparticle spectrum 

Following the method [12j used by us previously [3^, fS^jIB] we base the analysis 
of the Hamiltonian eq. ([1]) on getting back to the fermion (spinon) represen- 
tation by the standard move: 

Si = ^4^(y-o,pCip, (2) 

where cf^{cia) are the fermion creation (annihilation) operators; (Tafs are the 
elements Pauli matrices and the summation over repeating indices is as- 
sumed. Applying standard technique as described in Appendix [X] we reduce 
the problem to the set of 2 x 2 eigenvalue problems for each wave vector k: 

^k Ak \ / Uk \ _ p / Uk 
A^ -^k J[v^ J 'V ^k 
with 

Ck = -3 J^ Jt^t cos(kr) 

r 

Ak = 3 ^ JrAr cos(kr) 

T 

(summation over r extends to Tf, i = 1 -^ 4) which results in the excitation 
spectrum of the form: 



^k =± \/e^ + |A 



whose eigenvectors are combinations of the destruction and creation opera- 
tors with the above Bogohubov transformation coefficients. The above set of 
equations closes by the selfconsistency conditions of the form: 

e. = -^E-P(^kr)^tanh(§) 
A. = ^E^^P(-^kr)^tanh(^') (3) 

k 

for the order parameters (OPs) ^^, A^. The lattice symmetry considerations 
allow us to restrict ourselves by six OPs: ^a, ^c, ^ac] ^a, ^c, ^ac- Using the 
standard moves foreseen for the SU{2) symmetric solutions as described in 
Appendix \M we arrive to the quasiparticle spectrum: 

El = 9 {J^,C cos' x + j',CcCos^z + AJlCcCos'x cos' z), (4) 

where we set x = k^^; z = k^ and introduced effective OPs (a = a/Co + Va ' 

Cc = \fWT^c^ and Cac = \/ilc + vie ■ 

The spectrum eq. (jlj) is depicted in Fig. [2l 
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ID-RVB Cc 7^ 0; Ca, Cac = QID-RVB (c, Cac 7^ 0; Ca = 2D-RVB (c, Ca 7^ 0; C 







2D-RVBCc,Ca,Cac7^0 



Figure 2: Dispersion laws of the c-a-ca- RVB model for several exemplary values of the pseudogap/bandwidth 
parameters A = 3Ja(a] C = SJcCc] B = SJacCac indicating characteristic feature of the quasiparticle spectrum 
in different RVB states and the sketches of the relevant qDoS (see text for the details). 



If either of the OPs (a or (^ is the only nonvanishing OP, the quasipar- 
ticle spectrum acquires corresponding hues of nodes 2; = ±| or respectively 
X = ±^ where the quasiparticles have zero energy. Since the dispersion of 
quasiparticles takes place in only one crystallographic direction (a or c) we 
spell these states as one-dimensional (ID-RVB) states. In the low energy 
range the quasiparticle density of states (qDoS) in the ID-RVB states is con- 
stant. Due to the dispersionless ridge in the spectrum the qDoS diverges at 
the ceiling of the quasiparticle band [3]. This type of behavior is similar to 
that of the spinous of the anisotropic triangular Heisenberg model as found 
in [13] and considered in our previous work [31 lU |5l |6] in relation to CuNCN. 
If both OPs (a,c vanish and the OP (ac does not two paris of nodal lines exist 
along which the quasiparticles have zero energy. In this satate the qDoS 
diverges logarithmically at zero energy. We, however, do not focus on this 
pecular state and it is not shown in Fig. [2l 

If either of the nonvanishing OPs (a,c is complemented by the nonvan- 
ishing OP (ac quasi-lD-RVB (QID-RVB) states appear. The difference with 
the true ID-RVB states is that in the QID-RVB states there exists a nonvan- 
ishing dispersion in the direction transversal to the node lines. At the higher 
energies the quasiparticle spectrum of these states has maxima and saddle 
points instead of the ridge. Thus the qDoS develops a finite hop at the ceil- 
ing of the quasiparticle band and a van Hove singularty at somewhat lower 
energy, which serves as a pseudogap. By contrast, if the nonvanishing OP's 
are (c and (a then irrespective to the value of the OP (ac there are no lines 
of nodes, but four nodal points (k = (±^,±^)) of vanishing quasiparticle 
energies. In the vicinity of these points the quasiparticles are massless as in 
the disordered (pseudometallic) phase of the graphite monolayers [TJ]. The 
two possible states of this type are spelled as 2D-RVB ones. The qDoS in 2D- 
RVB states vanishes at the zero energy, being proportional to the energy well 
below the smaller pseudogap. Otherwise the quasiparticle dispersion law has 
saddle points at the two pseudogap energies and thus the qDoS of 2D-RVB 
state develops two van Hove singularities at the corrsponding pseudogaps. 



2.2 Free energy and phases of the model 

Following Ref. [10] one can write immeditaly the free energy of the c-a-ca 
model in terms of the above OPs: 



F = 3U: + 3 JcC + QJacC - ^ / In ( 2 cosh ( ^ ) ) d'k. (5) 



BZ 



where 6 = ksT (BZ stands for the integration over the Brillouin zone). Min- 
ima of eq. dS]) with respect to C's correspond to various possible phases of the 
system. We postpone the study of the ground state (zero temperature limit) 
of the present model to further publications and focus on the results which 
can be obtained with use of the high temperature expansion (technicalities 
are explained in Appendix [B|. 

The results obtained by minimizing the high temperature expansion of 
eq. ^ with respect to the OPs are depicted in the parameter phase diagram 
Fig. [3] and present in Table [1] Fig. [3] represents the phase diagram for the 
triple of exchange parameters subject to the condition Ja + Jc + Jac = 1 and a 
series of temperatures between 0.4 and 0.01 (the fractions of the above sum 
of the exchange parameters is meant). One can see that eight phases are 
possible. For the temperature above either of three critical ones: 

Of' = Pr. (6) 

the Curie paramagnetic phase: (1 - grey) in which all three OP's equal to 
zero persists. Below these temperatures gapless phases appear in respective 
corners of the parameter phase diagram. Their spectra are however quite 
different. There are two ID-RVB phases with OPs (cOT(a 7^ (2 - red 
and 3 - green areas in Fig. [3]). Due to the constant qDoS at the zero 
energy they have to exhibit a temperature independent paramagnetism which 
could be termed as Pauli paramagnetism had it been due to conductivity 
electrons. Thus we spell it as a quasi-Pauli one. The paramagnetic phase 
with (ac 7^ (4 - blue area in Fig. [3]) is quite different: the maximum of the 
spectrum at the ceiling of the quasiparticle band assures a constant value of 
the density of states (a finite hop) rather a divergence as in ID-RVB states, 
and logarithmically diverges at the zero energy. Thus one has to expect the 
paramagnetic susceptibility to logarithmically diverge at zero temperature in 
this phase. 
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Figure 3: Parameter phase diagrams for the c-a-ca-RVB model in the high- 
temperature approximation. The abscissa and ordinate represent respec- 
tively Ja and Jc, the entire parameters set is subject to the condition 
Ja + Jc + Jac = 1- The temperatures as fractions of Ja + Jc + Jac from the left 
upper diagram to the right lower one are 0.4, 0.35, 0.2, 0.15, 0.1, 0.04, 0.03, 
0.02, 0.01 (follow rows). The colour coding for phases is: Curie paramagnetic, 
all OPs are zero - gray, only one nonvanishing OP: red, green (ID-RVB, Pauli 
paramagnet), blue for CciCaiCac 7^ 0, respectively; one OP vanishing: mag- 
neta, cyan (QID-RVB, Pauli paramg-gnet), orange for (^q, otCo or^ac = 0, 
respectively - observe the order of the list; yellow codes the phase with three 
nonvanishing OPs. Orange and yellow phases (2D-RVB) feature combination 
of the gapped Arrhenius-like temperature dependence of magnetic suscepti- 
bility and its linear dependence well below the pseudogap. 



Below critical temperatures: 



ncrit '-^t' I i ^Jt' \ oJtJt' 



-1 



whose expression is already familiar from the high-temperature mean field 
analysis |3i i4j of the RVB states in the triangular anisotropic Heisenberg 
model, respective phases with two nonvanihsing OP's C,j. and C^^-i appear (the 
notation refers to transition from the state where only one OP C,- 7^ to a 
state where two OP's Cr? Cr' 7^ 0). The phases with two nonvanishing OPs 
as well are different. If one of the nonvanishing OPs is Qac (5 - magneta and 
6 - cyan areas in Fig. |3]) the qDoS is constant at the zero energy and thus 
quasi-Pauli paramagnetism has to be expected in these QID-RVB phases. 

The phase with three nonvanishing OP's (8 - yellow) is a transient one. 
It first appears below the octal point J a = Jac = </c = |; ^* = 1/^ where the 
Curie paramagnetic phase (grey) completely disappears and shows up from 
the above QID-RVB phases (magenta or cyan) at the critical temeperatures 

"M T- 

r,ac^a,c,ac " g ^^j^^ _ ^J^) ' ^ ^ 

(f = c, a for T = a,c) but exists above the critical temperature of: 



"i T T T 






where it switches to the 2D-RVB phase with only two nonvanishing OP's 
CcandCa (7 - orange). That latter phase appears also from the ID- RVB 
phases (red and green) at the critical temperatures given by eq. ([7]). The 
only difference between the dispersion laws of these two phases as depicted 
in 3-rd and 4-th columns of Fig. [2] is somewhat more pronounced disper- 
sion along a "ridge" in the case of the phase with three nonvanihsing OPs. 
Otherwise both 2D-RVB phases have a qDoS with two van-Hove singular- 
ities at the energies of their characteristic pseudogaps and the physics of 
these latter two phases has to be pretty similar. It can be checked that the 
transitions are largely of the second order, that is to say that the OPs split 
from zero continuously at the corresponding transition temperatures. The 
interphase borders between the ID-RVB and QID-RVB phases (2/5 and 3/6 
or red/magenta and gren/cyan, respectively) are special. Phases 2 and 3 
are unstable towards developing the nonvanishing OP (ac when the exchange 
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parameter Jac turns to be larger than, respectively, the exchange parameters 
Jc or J a at whatever temperature. On the lines Jac = Jc or Jac = J a where 
the OP Cac bounces from zero the transition temperatures from the Curie 
paramagnetic state (phase 1 - grey) to either ID-RVB or QID-RVB Pauli 
paramagnetic phases pairwisely coincide (that is to 2 and 5 i.e. to red and 
magenta or to 3 and 6 i.e. to gren and cyan). That means that if the system 
parameters fall in the corresponding (magenta or cyan) areas and it is cooled 
below the critical temperature eq. ([7]) with t = ac^r' = a,c , it directly goes 
from the Curie paramagnetic phase to the corresponding QID-RVB phase, 
rather to a ID-RVB phase. 

In the last column of Table [T] we show the analytical forms of the qDoS 
characteristic for the respective specific forms of the quasiparticle spectrum. 
It turned out quite unexpectedly, that these qDoS can be found analytically 
for all phases of the c-a-ca RVB-model. Leaving the details of the derivation 
for further publications we provide a sketch of the derivation in Appendix O 
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Table 1: Temeperature dependencies of the OPs for possible phases of the c-a-ca- RVB model in the high- 
temperature approximation. The areas of the existence of the corresponding phases are those where the 
expressions under the square roots are positive. The last column gives the qDoS in the respective phases. 
There K stands for the complete elliptic integral of the first kind. Details of derivation will be communicted 
elsewhere. Their characteristic graphs are given in Fig. [2l 
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3 Physical properties within the model 

3.1 Magnetic susceptibility 

We use the standard definition of tlie magnetic susceptibility per spin: 

,2 f^._.df{f 



X = -2/x^y g{e)-^de (10) 

wliere /(...) is tlie Fermi distribution function. Tlie qualitative behavior of 
the susceptibility as derived from the characteristic features of the qDoS 
is designated in the first column of Table [1] It is not surprizing that the 
temperature independent paramagnetism (conditionally denoted as "Pauli", 
although it goes about some other quasipartiles rather the band electrons in 
metals) takes place in the phases (Nos 2, 3, 5, and 6) with a constant density 
of states at the zero energy. The gapless phase 4 with two pairs of intersecting 
node lines manifests a logarithmic singularity in the qDoS at the Fermi level. 
This singularity, however, integrates and produces a logarithmic divergence 
of the susceptibility at the zero temperature. For two phases (Nos 7 and 8) 
with nodal points in the dispersion law and linear dependence of the qDoS 
in the low-energy range one has to expect as well a linear dependence of the 
susceptibility on temperature in the low-temperature region (well below the 
lower pseudogap) superimposed with a quasi-Arrhenius bahavior with char- 
acteristic energy of the pseudogap at higher temperature. That rich variety 
of possible phases on the c-a-ca-RVB model allows to eventually explain the 
magnetic behavior of CuNCN. 

As we mentioned above and previously [3], HI |5] the absence of magnetic 
scattering in CuNCN is perfectly explained by the hypothesis of the RVB 
character of its phases. The temperature independent paramagnetism of 
CuNCN as observed at higher temperatures is explained by formation of one 
of many ID- or QID-RVB phases (Nos 2, 3, 5, and 6). That means that 
only one of the OPs Q ot Qc is nonvanishing. As previously we assume that 
this phase sets on at some pretty high temperature which cannot be directly 
checked due to decomposition of the material P]. For a pseudogap to open 
in a (Q) ID- RVB state with, say, nonvanishing (^c (phase 2, red or phase 5, 
magenta) one needs that at an observable critical temperature the OP C,a 
splits from zero (phase 7, orange or phase 8, yellow) since for the pseudogap 
to open it cannot be the Cac OP. Thus we assume the following form for the 
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Figure 4: EPR magnetic susceptibility of CuNCN Ref. |1] (blue dots) as 
compared with the results of numerical integration with the qDoS for the 
phase 7. One can easily see the linear tail in the low temperatue range with 
no experimental points. Otherwise the parameters are 6^"^ = 95 K, z/ = 0.6, 
Ao = 145 K,C = 1250 K (red line) and r"* = 100 K, u = 0.5, Aq = 150 K, 
C = 1250 K (green line). 



temeperature dependence of the pseudogap A: 

e 



A{e)=AAl- 



Qcvit 



(11) 



below the critical temeperature ^"'* (where by ^"^^ one of the temperatures 



e 



crit 



or 9' 



icrit 
c,ac-^c,a,ac 



is meant) and perform the numerical integration of eq. 



(fTOj) for the susceptibility with the qDoS for the RVB phase with two pseu- 
dogaps where we also set S = Oo 

We performed several attempts and concluded that the values of ^™* 



{e\ 



icrit 



or 9 



crit 
c,ac^>c,a,acJ ^ 



Aq, C, and V can be adjusted to reproduce the exper- 



imental run of the susceptibility [4] . The results are shown in Fig. |H The 



^We notice that due to the character of the dependence of the OP Cac on the model 
parameters shown in Table[T]and the plausible assumption of the relation Jac ~ Jc between 
the exchange parameter, which needs to hold for the QID-RVB phase (magenta) to appear, 
this OP can never be large. Thus the characteristics of the system are basically not affected 
by the specific value of the Jac since its contribution is scaled down by the small value of 

sac- 
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value of C ir rather stable. Its scale is given by the magnitude of the Pauli 
paramagnetic susceptiblity at a higer temperature and can be used to fit the 
parameters of the original model eq. ([1]). The sets of parameters 6'™'* = 80 
K, z/ = 0.75, Aq = 240 K or 9^"' = 100 K, u = 0.5, Aq = 150 K, as well 
as z/ = 0.7, Aq = 170 K; u = 0.6, Aq = 145 K are equally good in terms of 
describing the susceptibility. The fact that the classical exponent ^ coming 
from our simple high-temperature treatment allows for an acceptable fit of 
experimental susceptibility is pretty remarkable. We assume at this point for 
the sake of simplicity that the phase with the finite qDoS at the zero energy is 
the ID-RVB one with the nonvanishing (c- Then using the zero temperature 
limiting value of this OP Cc = - [3] we get J^ = 1310 K. This nicely agrees 
with the original estimate of [T] {ca. 1000 K) and is more or less supported 
by other sources [8j. That is also what one can intuitively expect relying 
on the Goodenough-Kanamori rules [7jo Using it in the high-temperature 
estimate for the critical temperature eq. (Q we obtain for the highest critical 
temperature (that of the transition from the Curie paramagnetic to the ID- 
RVB state) the value 6^"^ = 490 K, which lies fairly above the decomposition 
temperature [1]. In case the QID-RVB phase is assumed to be responsible 
for the Pauli paramagnetism we take Jac = 1400 K as a plausible estimate. 
With use of it the highest critical temeperature 6"^^ = 435 K is as well high 
enough. In that sence these estimates are consistent. For the exchange pa- 
rameter Ja we notice that its values is pretty stable. Using eq. i^ for the 
temperature 6'^'^*^ ^ = 95 K of the pseudogap opening in the ID-RVB phase, 
we arrive to a reasonable value of Ja = 550 K. On the other hand consider- 
ing this critical value as the critical temperature 9'^'^^^^^ ^ ^^ given by eq. ([8]) 
yields Ja = 560 K which shows pretty good consistency between two options 
of defining the Pauli paramagnetic phase. 

3.2 Structure manifestations of transitions between RVB 
phases 

3.2.1 Theory 

Previously [5], [6] we were able to relate the RVB OPs and the lattice constants 
through "magnetostriction" - the linear coupling of the structure parameters 
with the exchange parameters. These moves absolutely apply in the present 



^Our previous estimate [J of 2300 K looks out to be somewhat exaggerated. 
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model with three exchange parameters and respective OPs. We assume as 
previously a linear relationship between the exchange parameters Jj. and 
geometry parameters p\: 

Jr = JrO + Yl Jr,xPX- (12) 

A 

Following [IB] we assume that zeroes of p's correspond to a hypothetical 
structure the CuNCN crystal would have provided the exchange interactions 
J-r are turned off. Deformation of this hypothetical structure requires the 
elastic energy 

fiX 

for each nearest neighbor Cu-Cu pair. The observed geometry of the crystal 
where the exchange interactions result in formation of one of the RVB phases 
corresponds to the minimum with respect to px of the free energy eq. ([5]) to 
which the elastric energy eq. (IT3|) is added and the exchange parameters are 
replaced according to eq. f lT2l) : 



^ + E^^AP, = 0. (14) 

Using the special form of the RVB free energy F (eq. (ISl)) we arrive to the 
explicit expression 

- the sought relation between the RVB OPs and their structure manifestations 
{At are numerical coefficients: 3 for t = a,c; 6 for r = ac) which further 
generalizes the famous bond-length-bond-order relation to the RVB states. 
Details are explained in Appendix [Dl 

3.2.2 Synchrotron measurements on CuNCN as explained by RVB 
phase transitions 

Now we can formulate what one could expect in the structural studies pro- 
vided CuNCN undergoes transitions between various RVB phases. The van- 
ishing OP's do not have any effect on the crystal structure. On the other 



hand, a transition accompanied by splitting from zero of some OP will be 
manifested in the structure changes as prescribed by eq. (1151) . The sign of 
the effect is controlled by that of the magnetostriction parameters J'^^. The 
temperature dependence of the structure described by the parameters px is 
thus that of the relevant combination of the squares of the corresponding 
OPs, which can be different in different phases. Previously E] we could 
relate the experimental data on the anomalous temperature dependence of 
the lattice parameter a with the ID-RVB to 2D-RVB phase transition of 
the anisoropic triangular Heisenberg model as accompanied by the opening 
of the gap in the quasiparticle spectrum in the a-direction. However, the 
experiment [H] showed some anomalies for the lattice parameter c at the 
temperature of the tentative ID-RVB to 2D-RVB phase transition (80 -t- 100 
K) and some more irregularities either in the a- or in the c-directions at ca. 
30 K. We also already mentioned above that the system of the exchange pa- 
rameters of the anisoropic triangular Heisenberg model was not particularly 
intuitive. 

The results of the structural studeies related to the lattice parameter 
a are shown in Fig. O We do not directly apply the general formula eq. 
( TTSj) . but, first, take into account the approximately diagonal form of the 
inverse matrix of the force constants as derived in Appendix [E] and thus 
conclude that the structure manifestations in the a- and c-directions are 
independent. Since the exchange parameters are sums of contributions of 
many superexchange paths one can expect that they depend on all types 
of interatomic separations that is the magnetostriction constants J'^^ with 
T = a,c, ac and A = a, c are not vanishing for all combinations of r and A. 
It seems, however, to be reasonable, that the exchange constant Ja is not 
dependent on the deformation in the c-direction and thus J'ac = ^- Under 
these assumptions the structural manifestations in the a-direction decouple 
and can be recovered by using r, \ = a, pa = Sa = (a — Qq) in eq. (ITSjl . 
Like previously [5l [6] we conclude that the equilibrium value of Sa in the 
RVB-phase with two pseudogaps must be proportional to the squared OP 
(a- Due to its direct relation with the pseudogap measured in the EPR 
experiements [3] we can relate two observed quantities: the deformations 
and the temperature dependent activation energy (the pseudogap ^4(6*)) in 
the quasi- Arrhenius regime: 

Sa = :^^A{9r (16) 
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Figure 5: The temperature dependence of the a-lattice constant as extracted 
(a) from synchrotron experiment [6j (red dashes with error bars) as confronted 
with the (^ values extracted from the ESR experiment and multiphed by a 
suitable constant J^ a/^aa = 0.08 A (blue dots) or confronted with the linear 
model eq. i^B) with the parameters^^™* = 100 K, u = 0.5, Aq = 140 K 
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(b) 



with the parameters 6™^ 
extracted from one of the magnetic susceptibility fits; (b) from the neutron 
experiment [17j (red crosses with error bars) and confronted with the linear 
model eq. (IT6|) . 



The crucial point is the sign of the effect. At the first glance the situation 
seems to be counterintuitive since in order to conform with the experiment 
(increase of the a-parameter in the 2D-RVB phase with two two-pseudogaps 
as compared to the (Q)ID-RVB phase (2, red - or 5, magenta) the exchange 
parameter Ja has to increase with the increase of the Cu-Cu interatomic 
separation. However, the effective value of the (antiferromagnetic) exchange 
parameter Ja is a sum of numerous contributions of different signs: 

Ja = t/a(antiferro) — Ja(ferro) > 0, 

where both Ja(antiferro) and Ja (ferro) are positive. The antiferromagnetic 
contribution is accumulated by summing up contributions from numerous 
superexchange paths and is weakly affected by the Cu-Cu interatomic sepa- 
ration since no direct matrix elements between the states of the two Cu atoms 
affects Ja(antiferro). By contrast the ferromagnetic contribution strongly de- 
pends on the angle CuNCu (see Refs. [TJ EI]) and in the range of CuNCu > 
90°, which is the case for CuNCN, decreases while the CuNCu angle and thus 
the Cu-Cu separation increases. Respectively the effective antiferromagnetic 
exchange parameter Ja increases while the counterpoising ferromagnetic con- 
tribution decreases. This explains the overall positive sign of J^ ^ and thus 
the experimental fact of increase of the lattice parameter a in the phase with 
the pseudogap developing in the a-direction. 

The quantitative agreement (shown in Fig. [5](a) by blue and green dots) 
is achieved by confronting the amplitude of the structure effect 4 ■ 10^"^ A [6] 
in the a-direction with the zero temperature limit of the pseudogap A{9 — t- 0) 
of ca. 70 K j3]. The problem with the latter estimate is that it is determined 
within somewhat different model [U [5]. Alternatively we take the value of 
Aq = 140 K and the corresponding classical value of the exponent z/ = | as 
a plausible estimate for the zero temperature limit for the pseudogap and 
with the above value Ja = 560 K we get (a{(^ — > 0) = 0.08. Then the 
ratio of effective constants finally responsible for the spin-phonon coupling is 
J'aa/^aa = 0.19 A. Combining this with the value of the force constant Kaa 
derived in Appendix [El we can now evaluate J^^ to be ca. 12000 K/A which 
seems to fairly fit the expectations. 

The synchrotron experiment was subject to a heavy and unfair criticism. 
Thus the neutron difraction study of the temperature dependence of the 
CuNCN structure have been performed [T7]. The results concerning the a 
lattice parameter are shown in Fig. El^b). The most remarkable is that the 
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amplitude of the structure change as measured in the neutron scattering co- 
incides with that coming from the synchrotron one: 4 ■ 10~^ A. That means 
that our previous estimates of the parameters rationahzing the synchrotron 
and ESR data remain vahd. In Fig. [5t^a), (b) the experimental data fairly 
lie on a straight line in the temperature range between 100 K (the mea- 
sured minimum of a determined by the neutron difraction and somewhat 
less characteristic for the synchrotron data) and 30 K. This agrees with an 
assumption that in a rather wide range the temperature dependence of the 
OPs follow the standard temperature behavior with the classical value of 
the critical exponent i' = \ since accoridng to our treatment eqs. (IT5|) . (IT6|) 
the lattice parameters must obey the temperature evolution with the expo- 
nent 2z/ i.e. be linear, as observed. This, however, changes to somewhat 
more chaotic behavior at ca. 30 K. This incidentally corroborates with the 
temperature run of the c lattice parameter measured in the synchrotron ex- 
periment ^ (Fig. [6]) which did not receive due attention so far. As one can 
see the lattice parameter c generally (almost linearly) decreases with temper- 
ature. It shows some irregularities (not changing the sense of the run) about 
80 -i- 100 K, where the magnetic susceptibilty changes from the quasi-Pauli to 
the quasi- Arrhenius regime due to opening the pseudogap in the a-direction. 
However, the sense of the temperature run of c changes at 30 K where the 
lattice parameter a starts to show irregularities as well. 

In order to rationalize the low-temperature behavior of the lattice pa- 
rameter c and eventually of a we introduce one more geometry variable 
2pc = 5c = c — Co and apply the general relation eq. (fT5|) and by this show 
that the variation of two lattice parameters is: 

6C J \K\ \ {-K,cJ'a,a + 2Ka,aJ;.) C + [-KacJ'c,a + 2i^a,aJ;.) Cc 

(17) 

where we denote by \K\ the determinant of the 2x2 matrix of the force 
constants and make use of our previous consideration allowing to neglect the 

OP Cac 

In order to qualitatively understand the temperature run of the lattice 
parameters a and c formally given by eq. (IT7I) we take a closer look at the 
lower-temperature range of the parameter phase diagram Fig. [3l There one 
sees that below the octal point the thermal evolution consists en gros in 
squeezing out all the phases by the 2D-RVB phase with vanishing Cac (7 - 
orange). The 2D-RVB phase with three nonvanishing OP's (8 - yellow) is 
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Figure 6: The temperature dependence of the c-lattice constant as extracted 
from the synchrotron experiment [6]. 
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Table 2: Temeperature dependencies of the OPs for possible phases of 
CuNCN in the high-temperature approximation of the c- a- ca- RVB model. 
The areas of the existence of the corresponding phases are given by the con- 
dition of positiveness of the expressions under the square roots. 



Color code 



OP's vs. 9 



Magnet a 





3Jc 



1- 



cnt 
c,ac 



9. 



oJacJ c 



c,ac 8{3Jac-2Jc) 



2 e e{Jac-Jc) 

"J •Jac V •Jac'Jc 



Yellow 



Ca 
Cc 
Sac 



^.-n 



3Jc 



icrit 
c,ac—>a,c,ac 



^crit __ 3JacJa 

"c,ac^a,c,ac S{3Jac-2Ja) 



1 6> 

3 Jac 



cnt 
a,aC'^a,c,ac 



cnt 
a,c,ac 



9 



crit 



'J'J acJ c 



a,ac^a,c,ac 8{3Jac-2Jc) 

1 C*ac = eq.(l9D 



Orange 



461 



15J„ 
46> 
VTSJc 

= 



cnt 
c^fa,c 



/icrit __ 3JaJc 
^c^a,c 8{3Jc-2Ja) 



crit 



9' 



OjgJc 



a^a,c "a^a,c 8(3Ja-2Jc) 



transient, its parameters' area is never large and it is subject to deformations 
and displacements under the "pressure" of the phase 7. The same "pressure" 
squeezes the area of existence of the Q ID- RVB phase 5 (magenta) as well. 
Assuming the position of CuNCN on the parameters phase diagram Fig. 
[3] within the magenta phase, but close and somewhat below the quadruple 
point of the phases 2, 5, 7, and 8 (red, magenta, orange, and yellow) in the 
rightmost graph in the middle row one can see on the following slices the 
following sequence of transitions between the RVB phases: 



QID - RVB 


-)■ 


2D - RVB 


-)■ 


2D - RVB 


S>c; Sac 


— )■ 


(,c; S>ac) Sa 


-)■ 


CcCa 


magenta 


— )■ 


yellow 


-)■ 


orange 



The thermal dependence of the OPs within these phases is described by the 
formulae of Table |2l One should not expect that these expressions derived 
from the high-temperature expansion for the free energy are exactly valid at 
low-temperature. Specifically, the prefactors 9 should not be taken seriously 
since we expect that the C^a,c OPs flow to some finite values as temperature 
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flows to zero, although splitting from zero vaues at the critical temperatures 
given next to them. Upto that uncertain factor the temperature dependence 
of the OPs in three presumably observed phases is as follows. As we explained 
earlier, since the exchange parameters satisfy the condition Jac ~ Jc the OP 
C,ac is always small and decreases with decreasing temperature, thus we do not 
consider it explicitly further. The OP Cc should in principle increase either 
in the magneta or in yellow phases, but we assume that in the interesting 
temperature range 9 ^ 6'™* = 6'™*^_j^^ cac — 435 K it has almost reached its 
zero temperature limit and does not significntly change any moreo The OP 
C^a splits from zero at the critical temperature of 0'^^^^^^^ cac^ 100 K (transi- 
tion to the yellow phase) which affects the a lattice parameter as explained 
above. In the general setting eq. f ll7p . but under previous assumption of 
J'ac — ^ the temeperature independent (^ does not contribute to the visible 
variation of the lattice constants. However, the lattice parameter c turns 
out to be sensitive to the pseudogap opening in the a direction through the 
off-diagonal element of the the inverse matrix of the force constants. Partic- 
ularly remarkable is the fact that the effect on the c is predicted to have the 
sign opposite to that on a (a increases, c decreases), as observed. Finally at 
the lowest accessible critical temperature ^^"^ac ^q- & a transition to the or- 
ange phase takes place. A word of caution needs to be said here: the present 
set of the exchange parameters yields very narrow temeperature range where 
the transient (yellow) phase can exist. We assume that eq. (^ strongly over- 
estimates this temperature (92 K), which must by considerably lower {ca. 
30 K) since one cannot rely upon the results of the high-temperature ex- 
pansion any more. For obtaining more reliable estimates one needs to know 
at least the zero temperature limit of the a-c-ca model which is yet to be 
done. However, we assume that the general shape of the phase diagram is 
correctly reproduced by the high-temperature expansion and proceed within 
this setting. The evanscence of (ac at 0^^^^^ whatever it is affects neither the 
bandwidth, since this OP is never large, nor the character of the temperature 
dependence of (a, although the variation of the slope can be as well expected. 
However, at this phase transition a remarkable change can be expected to 
the character of the temperature dependence of (c- Namely, it switches from 
increase to decrease. This happens because of an instanteneous change of 



■^Possible minor effect of increase of Cc produces no visible temperature dependence of 
the QID-RVB quasiparticle bandwidth due to decrease of (ac since these two contributions 
changing in opposite directions compensate variations of each other. 
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crit 
a—>a,c 



its reference temperature shown in Table [2] from positive ^™ c-s>a c ac ^° ^\ 
which is negative (and much smaller by absolute value: ca. —300 K) in the 
relevant area of the exchange parameters' {Jac ^ Jc > Jo) space. Together 
changing the the sign and the magnitide ot the reference temperature causes 
the change of the temperature run of the OP Cc which starts to decrease in 
the orange phase and through eq. fll7p affect the temperature behavior of the 
lattice parameters a and c. This tentatively and qualitatively explains the 
anomalies observed in Figs. [5] and [6] although the available amount of experi- 
mental data does not suffice to reliably estimate quite a few magnetostriction 
and other parameters required for at least semiquantitative description. In 
the minimal setting: i.e. neglecting the off-diagonal elements of the inverse 
matrix of force constants and the off-diagonal magnetostriction terms we 
immeditely obtain 

which allows to at least conclude that the magnetostriction parameter J^^ 
has intuitively understandable sign: it is negative since 5c turns to be positive 
when the variation of OP Cc is negative. 

4 Conclusion 

A new form of the frustrated spatially anisotropic antiferromagnetic Heisen- 
berg Hamiltonian close to the popular J1J2J3 model with exchange param- 
eters Jc, Ja, and Jqc extended along the c, a, and a±c directions of a two- 
dimensional rectangular lattice is proposed. When applied to model fas- 
cinating physics of copper carbodiimide (CuNCN) it explains the absence 
of magnetic order in CuNCN down to 4 K by assuming resonanting va- 
lence bond (RVB) character of the emerging phases. The quasiparticle spec- 
trum of the RVB model of the proposed Hamiltonian shows three princi- 
pal regimes: (i) a state with two pairs of lines of nodes, (ii) states with 
a pair of lines of nodes (termed as ID- and QID-RVB states), (iii) states 
with two pseudogaps and four nodal points (2D- RVB states). Extraordi- 
nary rich parameters-temperature phase diagram of the model contains eight 
different phases whose magnetic behavior includes Curie and quasi-Pauli 
paramagnetism (ID- and QID-RVB phases), and (pseudo)gapped (quasi- 
Arrhenius) paramagnetism (2D-RVB phases). Adding magnetostriction and 
elastic terms to the free energy of the model explains the the temperature 
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dependence of the CuNCN crystal structure by assuming that a sequence 
of transitions between different RVB phases occurs in CuNCN while tem- 
perature decreases. Confronting the model with the magnetic susceptibility 
and strucutre (both synchrotron and neutron) data recorded as functions of 
temperature in the range between ca. 10 and 200 K shows a remarkably good 
agreement between the theoretical predictions and the experiment which is 
reached by ascribing the values of the model parameters which are intuitively 
acceptable both in terms of their absolute magnitudes, relative values, and 
the character of their geometry dependence. 
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A Equations of motion and self consistency equa- 
tions. 

Equations of motion are based on the Heisenberg representation in which 
each operator obeys the following: 

ihA = [A, H] (18) 

where [...,...] stands for the commutator and "'" for the time derivative. Ap- 
plying this to the creation and annihilation operators c^^{cra) and performing 
commutation, mean field decoupling and Fourier transformation as done pre- 
viously |3] results in mean field equations of motion for them: 

3 V 

T 

3 V 

^^Ci^a = 2 5Z ^^^^ ^°^(^^)^k/3 -fEr^^^T cos(kr)c_k_^ (19) 

r 

These latter reduces to the set of 2 x 2 eigenvalue problems for each wave 
vector k: 



^k ^k \ f ui,\ _ f Uk 
K -^k j[vkj ' Uk 



with 



Ck = -3 ^ ^T^T cos(kr) 

T 

Ak = 3 ^ JrAr cos(kr) 

T 

(summation over r extends to Tj; i = 1 ^ 4) which results in the eigenvalues: 



Ek =± V^k+|A 

whose eigenvectors are combinations of the annihilation and creation opera- 
tors with the coefficients Wk, f k- This set of equations closes by the selfcon- 
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sistency conditions: 

e. = -^E-P(^kr)^tanh(^) 

k ' 

^^ = 45:-p(-*-)t*^""(t) '"" 

k ^ ' 

for the order parameters (OPs). The lattice symmetry considerations alfow 
us to restrict ourselves by the OPs: C,a,^c,^ac]^a,^c,^ac- For the complex 
OPs we introduce a polar representation: 

The standard moves foreseen for the SU{2) symmetric solutions are used 
to exclude the cross terms in OPs from E^ which leads to the system of 
equations: 

(three equations). Introducing the relative phases as: 6a = fa — fac] dc = 
fc ~ fac we arrive to the equations of the fornio 

^a^c + VaVc COs{6a - 0^) = 
^aCac + VaVac COS 9a = 
iciac + rjcTjac COS 6'c = 

similar to those derived in [9] which can be satisfied e.g. by setting 

(21) 
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^0 
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lac 
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Of course, the OPs can be also vanishing. 

Under the above conditions the spectrum of quasiparticles acquires the 
form: 

K = 9 ( J2 {ea + vl) cos^ X + J2 {e + vl) cos^ z + AJl [il + r,l) cos^ x cos\2) 

where we set x = k^; ; 2 = k^ . 



^The consequences of setting Oc = fac — fc and thus having the first equation in the 
form £_a£.c + VaVc cos{9a + Oc) = 0, which exactly coincides with [9], yet have to be studied. 
Most probably it brings up more degenerate phases with various phase angles. 
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B High-temperature expansion. 

At high temperature we can use an expansion: 

which when inserted in eq. ([5]) integrates exphcitly. For determining the 
critical temperatures to the first approximation it suffies to restrict ourselves 
by the second power terms. This results in an expression quadratic in (t-. 
Combining thus obtained "kinetic" energy: 

J '1 -2^2 I 1 t2a2 , t2 a2 



with the potential energy terms from eq. ([5]) we get: 

9 
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j: + 3x 



Cr- 



This result can be improved with use of the Ginzburg-Landau approximate 
free energy -Fgl(Ct)^) which involves the terms up to the fourth power in 
(^T-'s. They appear from the integration of the 4-th power of the spectrum 
which performs explicitly and yields the "kinetic" energy of the form: 

^ [36CJ'. {27CJI + %lJl - 2Ae^) 

+ 486^71 - 2mHlJl + ^KUt] 

which together with the potential energy yields the free energy Fgl(Ct;6') 
used for further analysis. 

C Quasiparticle densities of states in various 
RVB phases 

In Section [2l2] we gave an impression of the complexity of the phase diagram 
of the RVB model with three exchange parameters. We also gave a brief 
description of the most characteristic features of the qDoS in various RVB 

32 



phases. Here we provide a brief sketch of the derivation of qDoS given in 
Table [H 



The definition of the qDoS reads: 

9{e) = ^l6{e-E^)d'k 



BZ 



Following Ref. [15.j we insert the integral representation for the Dirac 6- 
f unction: 

oo 

9{e) = ^^ J J dtexp {it {e - E^)) d^k = 



BZ -oo 

oo 



1 . , 

. dt 

271 



exp iite) — - / exp i—itE]^ d k 
A-n^ J 

BZ 



For all phases having lines of nodes (those with numbers 2-6) the integration 
over one of the components of the wave vector k in the BZ is performed and 
yields an intermediate result in terms of the the Bessel and Struve functions 
of arguments dependent on the Fourier transformation variable t and the 
remaining component of the wave vector k. The Fourier trasnforms with 
respect to t can be done for the intermediate answers of that form. It yields 
integrands of the elliptic integrals over the remaning component of the wave 
vector k. This solves the problem of calculationg the qDoS for the dispersion 
laws with the lines of nodes in the BZ. The results are given in respective 
cells of Table [1] 

In order to cope with remaining two phases whose spectra contain only 
nodal points (numbers 7 and 8) one more trick, namely performing previ- 
ous moves for the squared spectrum and thus obtaining the distribution of 
the states as a function of their squared energy solves the problem. The 
distribution of squares of the quasiparticle energies is: 



si^')-^.ji{^'-El)<e^. 



BZ 

Then the sought qDoS is given by |16| : 

g{e) = 2eQ{e^) 
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To obtain ^(e^) we again use the integral representation of the Dirac S- 
function: 



g{e 



1 1 



oo 

/ / rftexp {it {e^ - El)) cfk 



BZ -oo 



27r 



dt 



exp {ite^) 



A-K^ 



exp {-itE"^ <f\i 



BZ 



Remarkably enough sequential integrations of the squared spectrum over 
one of the components of the wave vector k and Fourier transformation with 
respect to t yield the expressions of the same form as integration of the 
spectrum itself: the Bessel function and an integrand of the elliptic integral. 
Thus the final intergation over the remainig component of the wave vector k 
yields some elliptic integrals given in respective cells of Table [1] 

D Theory of the structural manifestations of 
the RVB states. 

We start from the mechanic equilibrium condition of the crystal in a RVB 
phase: 

Y.K^^p, = Q. (23) 



dp> 



fiX 



and notice that the free energy F given by eq. ([5]) has a special form. The 
first term ("kinetic" energy) is integral of a function of the dispersion law £'k 
whose argument has the form: 

T 

where Ar are numerical coefficients (3 for r = a, c; 6 for r = ac); fr (k) are 
trigonometrical expressions (cos k^ for t = a,c; cos k^ cos kc for r = ac). The 
"potential" energy contribution to the free energy eq. (^ is: 



34 



Due to the above special form of the "kinetic" and "potential" energies the 
self-consistency equations for (j. have the form: 

I tanh (^^^ E'^AlJXfr (k) d'k = A^Ur 

BZ 

each, to be solved simultaneously for all r. Apparently whatever (sub)set of 
Cr = satisfies the equations. For the nonvanishing OPs the self-consistency 
equations acquire the form: 



tanh^^jK/.^(k).'k=— . (24) 

BZ 

Now we can turn to the structure determination. The derivative of the 
kinetic energy with respect to the geometry parameters px reads: 

-2 Y. Al (j.CV;, + JX^j^ I tanh (^) E^ (k) d'k. 

^ BZ 

The value of the above integral for the equilibrium values of the OPs is given 
by the self-consistency conditions eq. (12^ . Thus the derivative of the kinetic 
energy rewrites: 

Combining this with the derivative of the potential enenrgy and that of the 
elastic energy we get: 

-2 5^ A^ (CJr,X + JrCr^) +E A^ U^X + ^JrCr^) +^1 ^ ,XP , = 

The terms including the derivative ^ stemming from the kinetic and poten- 
tial energies cancel each other: a remarkable consequence of the Hellmann- 
Feynman and virial theorems, which immediately results in: 

T flX 

which after some trivial algebra results in: 
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P. = EiKy.l[Y.^rCJr,A- (25) 

E Estimate of the force matrix from the elastic 
constants 

Now we notice that the deformation tensor (for methods used for this and 
further evaluates see Ref. [19j) corresponding to the structure variation as 
desribed by the geometry deformation parameter pa = Sa has only one non- 
vanishing component Uaa = 5a/ a^. The elastic energy of the unit cell under 
such deformation is: 



-Caa.aaflO^OCo — = 4 X -Kaa {^O) , 

2 V^o / 2 

where Caa,aa is the corresponding element of the elasticity moduli tensor, 
the multiplier of 4 in the right hand side takes into account that each unit 
cell of CuNCN contains four Cu-Cu interactions along the a-direction, and 
Oq, feo; and Co are the orthorombic (?) lattice constants, so that we get the 
estimate: 

K -b^c 

4ao 
Completely analogous consideration of the deformation parameter 2pc = (^c, 
yields the defomation tensor with single nonvanishing component Ucc = 5c/ Cq 
(the factor 2 takes care of the fact that in the structure shown in Fig. [T]the 
Cu-Cu distance in the c direction fits twice in the unit cell) yields: 

^r< „ I, „ (^(^\^ 1^ „ I, „ C^PcV ,1^ ao6o o a^At^ „2 



-CccccflO^OCo — = oC'ccccOO^OCo = 4X-Ccc,cc Pc = 4:X-KccPc 

^ \Co / ^ \ Co / ^ Co I 

where the factor 4 in the rightmost expresion as in the case of the lattice 
direction a takes into account the presence of the four Cu-Cu bonds in the 
ah cross section. Thus: 

K -^C 

Co 

In the general case when both lattice parameters a and c change although 
not b and neither of angles the deformation tensor has two nonvanishing 
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components Uaa and Ucc defined above. The energy of one unit cell under 
such deformation reads 



2 



, , 1^ f^A 1^ f^cy ^ f6a\ (5c 
2 V«o/ 2 \coJ \aoJ Vco. 

Singling out the remaining off-diagonal term we write: 

boCaa,ccSaSc = boCaa,cc'2paPc = 4 X 4 X KacpaPc 

where two factors of 4 in the rightmost term take into account the fact the 
four Cu-Cu bonds are extended from the given unit cell in either direction a 
or c. Thus 

K -—C 

o 

From the VASP calculations on CuNCN in various antiferromagnetic 
states we have for Caa,aa the estimates of 195.6 or 181.3 GPa, Ccc,cc is 159.6 
GPa, and Caa,cc 49.5 GPa [20j. Taking into account the SI units relations 
(1 GPa = 10^' J/m^, lA=10-^° m, ks = I.SS-IO^^^ J/K and the values of 
the lattice constants of CuNCN (ao = 2.99, bo = 6.19, cq = 9.41 A, see Refs. 
[H E]) we obtain the elastic constant Kaa to be in the range of 64100 to 
69000 K/A^ (i.e. ca. 6 to 7 eV/A^) pretty smaller than the characteristic 
values derived for analogous constant in polyenes ^E\, which one, however, 
could expect provided the difference between the concerned deformations of 
intramolecular bonds in polyene and somewhat weaker interionic interactions 
in CuNCN. Two other force constans get estimates Kcc = 22750 K/A^ and 
Kac = 2775 K/A^. With that large difference between the diagonal and 
off-diagonal force constants we can sometimes assume the matrix K~^ to be 
diagonal with elements equal to inverse diagonal force constants K~^ and 
K: 
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